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§1. Introduction 



The problem of predicting from first principles total cross sections at high energy is one of 
the oldest open problems of hadronic physics (see, for example, RefP and references therein). 
As QCD is believed to be the fundamental theory of strong interactions, it should predict the 
correct asymptotic behaviour: nevertheless, a satisfactory explanation is still lacking. The 
problem of total cross sections is part of the more general problem of high-energy elastic 
scattering at low transferred momentum, the so-called soft high-energy scattering. As soft 
high-energy processes possess two different energy scales, the total center-of-mass energy 
squared s and the transferred momentum squared t, smaller than the typical energy scale 
of strong interactions (|t| < 1 GeV 2 <C s), we cannot fully rely on perturbation theory. A 
genuine nonperturbative approach in the framework of QCD has been proposed by Nacht- 
mann in RefP and further developed in Refs.P^ using a functional integral approach, 
high-energy hadron-hadron elastic scattering amplitudes are shown to be governed by the 
correlation function of certain Wilson loops defined in Minkowski space. Moreover, as it 
has been shown in Refs.P - ^ such a correlation function can be reconstructed by analytic 
continuation from its Euclidean counterpart, i.e., the correlation function of two Euclidean 
Wilson loops, that can be calculated using the nonperturbative methods of Euclidean Field 
Theory. 

In RefsP^ 1 '^ we have investigated this problem by means of numerical simulations in 
Lattice Gauge Theory (LGT). Although we cannot obtain an analytic expression in this way, 
nevertheless this is a first-principle approach that provides (inside the errors) the true QCD 
expectation for the relevant correlation function. In this contribution, after a quick survey 
of the nonperturbative approach to soft high-energy scattering in the case of meson-meson 
elastic scattering, we will present our numerical approach based on LGT, and we will show 
how the numerical results can be compared to the existing analytic models. 

§2. High— energy meson— meson elastic scattering amplitudes and Wilson— loop 

correlators 

We sketch here the nonperturbative approach to soft high-energy scattering (see RefP 
for a more detailed presentation). The elastic scattering amplitudes of two mesons (taken 
for simplicity with the same mass m) in the soft high-energy regime can be reconstructed, 
after folding with the appropriate wave functions, from the scattering amplitude M.(dd) of 
two dipoles of fixed transverse sizes -Ru_, R2±, and fixed longitudinal-momentum fractions 
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/i, f 2 of the two quarks in the two dipoles: 31 

M (M) (s,t;l,2) = -i2s [ d 2 z ± e tq ^C M ( X ~ log( S /m 2 ); 1, 2), 



(2-1) 



where s = (pi + P2) 2 and t = — \q±\ 2 (q±_ being the transferred momentum) are the usual 
Mandelstam variables, and the arguments "1" and "2" stand for u R\x, fx" and u R2±,fz' 
respectively. The correlation function Cm is defined as the limit Cm = lim Qm of the 
correlation function of two loops of finite length 2T, 



Gm{x;T;z ± ; 1,2) = 



where (. . .) are averages in the sense of the QCD functional integral, and 

,(T) 



(2-2) 



-ig (p A ll {x)dx tJ ' 

'Ci, 2 



(2-3) 



are Wilson loops in the fundamental representation of SU(N C = 3); the paths are made up 
of the classical trajectories of quarks and antiquarks, 



d : X 1 ^ ( t ) = z + ^t + ff l R u C 2 : X 2 ^ (r) = + ff l R 



Mi] 



P2 



?i[q] 



(2-4) 



m m 
with t G [— T, T], and closed by straight-line paths in the transverse plane at r = ±T in 
order to ensure gauge invariance. Here 



Pi = m ( cosh ^ , sinh ^ , 



p 2 = m ( cosh — , — sinh — , Oj 



(2-5) 



X being the hyperbolic angle formed by the two trajectories, i.e., Pi ■ p 2 — m2 coshx. More- 
over, R 1 = (0, 0, R lA _), R 2 = (0, 0, R 2± ), z = (0, 0, z ± ), and ft = 1 - f u /? = (i = 1, 2), 
with /j G [0, 1] the longitudinal-momentum fraction of quark "z" . 
The Euclidean counterpart of Eq. (12-21) is 



g E (9;T;z ± ;l,2) 



{^ ( P)e{w { P) 1 



(2-6) 



where now {■ ■ )e is the average in the sense of the Euclidean QCD functional integral, and 
the Euclidean Wilson loops 



j/P exp 


-ig <b 




JCx 



are calculated on the following straight-line paths, 



d : .V;-' (r) = z + + frRiE, C 2 : X™ (r) 
m 



Pie 



Mq] 



^ + J 2 

m 



(2-7) 



(2-8) 
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with r G [-T, T] , and closed by straight-line paths in the transverse plane at r = ±T. The 
four-vectors pi E and p 2 E are chosen to be (taking X E4 to be the "Euclidean time" ) 

/ 9 9\ / 9 9\ 

Pie = m^sin-,0_L,cos-J, p 2E = - sin -, 0j_, cos -J , (2-9) 

9 being the angle formed by the two trajectories, i.e., Pie • P2E — m 2 cos9. Moreover, 
Rie = (0, R\±,0), R 2 e = (0,i?2±,0) and ze = (0, z±,0) (the transverse vectors are taken 
to be equal in the two cases). Again, we define the correlation function with the IR cutoff 
removed as Ce = lim Q E . 

T->oo 

It has been shown that the correlation functions in the two theories are connected by the 
analytic-continuation relations®~^ 

g M { X ; T- z ± ; 1, 2) = Q E {~ix\ iT\ z ± ; 1, 2), V X G M + , 
g E {9;T;zr, 1,2) = g M {i9;-iT;z ± ; 1,2), V^g(0,vt). (2-10) 

Here we denote with an overbar the analytic extensions of the Euclidean and Minkowskian 
correlation functions, starting from the real intervals (0, 7r) and (0, 00) of the respective 
angular variables, with positive real T in both cases, into domains of the complex variables 
9 (resp. x) an d T in a two-dimensional complex space. (See RefP for a more detailed 
discussion: in particular, in RefP we have shown, on nonperturbative grounds, that the 
required analyticity hypotheses are indeed satisfied, thus obtaining a real nonperturbative 
foundation of Eqs. (IgTOjl .) 

Under certain analyticity hypotheses in the T variable, the following relations are obtained 
for the correlation functions with the IR cutoff T removedP'^ 

C M (x; zr, 1, 2) = C E (-tx; z±\ h 2), V x g R + , 

C e (9; zi; 1,2) = C M (i9; z ± ; 1, 2), W G (0, tt). (2-11) 

Finally, we recall the so-called crossing-symmetry relaUons^-^ 

C M {m - X ; z ± ; 1, 2) = C M {x\ z±; 1,2) = C M ( X ; z x] T, 2), V X G M + , 

C E (ir - 9; z ± ; 1, 2) =C E {9;z ± ; 1,2) = C E (9; z±;l,2), V^g(0,vt). (2-12) 

Here the arguments V stand for U -R i± , 1 - f" (i = 1, 2): the exchange "1, 2" ->• "1, 2", or 
"1, 2" — y "1, 2", corresponds to the exchange from a loop-loop correlator to a loop-antiloop 
correlator, where an antiloop is obtained from a given loop by exchanging the quark and the 
antiquark trajectories. 

In the following, we will take for simplicity the longitudinal-momentum fractions f\, f 2 of 
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the two quarks in the two dipoles to be fixed to 1/2: as it is explained in the Appendix of 
Ref.p^ one can always reduce to this case without loss of generality. We will also adopt the 
notation Q E {9\ T; z x ; R\_l, R 2 ±) = G E (9;T; z±; R 1± , f x = ^,R 2 ±,f 2 = §), and similarly for 
C E . 

§3. Wilson— loop correlators on the lattice 

The gauge-invariant Wilson-loop correlation function Q E is a natural candidate for a 
lattice computation, but some care has to be taken due to the explicit breaking of 0(4) 
invariance on a lattice. As straight lines on a lattice can be either parallel or orthogonal, we 
are forced to use off-axis Wilson loops to cover a significantly large set of angles.^ To stay 
as close as possible to the continuum case, the loop sides are evaluated on the lattice paths 
that minimise the distance from the true, continuum paths: this can be easily accomplished 
making use of the well-known Bresenham algorithm 1 '?' 1 to find the required "minimal distance 
paths' 1 corresponding to the sides of the loops. The relevant Wilson loops Wl(/||; r±; n) are 
then characterised by the position n of their center and by two two-dimensional vectors In 
and r_j_, corresponding respectively to the longitudinal and transverse sides of the loop. 

On the lattice we then define the correlator 

yL{hhh\\;d±;ri ± ,r 2 ±) = ~ ^ NV ~ -p - 1, (3-1) 

(W L (/ni;f 1± ;rf))(>V L (/ 2 ||;f 2± ;0)) 

— * 

where d = (0, d±, 0), and, moreover, 

CL{lihl2\\',d±;fi±,f 2 ±) = lim GL(h\\,h\\ld±-,ri±,r 2 i_), (3-2) 

Li,L,2—>oo 

where Li = are defined to be the lengths of the longitudinal sides of the loops in lattice 
units, and Li = h\\/Li. In the continuum limit, where 0(4) invariance is restored, we expect 

GL(h\\,h\\;d±;fi_ L ,f 2 ±) ~ Qe{0]T x = aL 1 /2,T 2 = aL 2 /2;ad_ L ;ar 1 _ L ,af 2 ±), 
C L (L\,L\;d±-,r 1± ,r 2A _) ~ C E {0;ad±;afi±,af 2 i.), (3-3) 

where ■ l 2 \\ = cos 9 defines the relative angle 9 and a is the lattice spacing. 

To keep the corrections due to 0(4) invariance breaking as small as possible, we have kept 
one of the two loops on-axis and we have only tilted the other one as shown in Fig. HJ the 
on-axis loop Wli is taken to be parallel to the x E \ axis, lu = (L x , 0), and of length L\ = 6, 8. 
We have used two sets of off-axis loops Wl2 tilted at cot 9 = 2, 1, 1/2, 0, —1/2, —1, —2, i.e., 
9 ~ 26.565°, 45°, 63.435°, 90°, 116.565°, 135°, 153.435°. We have used loops with transverse 
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Fig. 1. The relevant Wilson-loop configura- 
tion. Using the 0(4) invariance of the Eu- 
clidean theory we have put p±E parallel to 
the xei axis. 



f =Mi + 1 } 



Fig. 2. Loop configuration in the transverse 
plane. In the li ave" case the link orienta- 
tion is not shown as it is averaged over. 



size |rxj_| = |r* 2 j_| = 1 in lattice units; the loop configurations in the transverse plane are 
those illustrated in Fig. [21 namely d±_ \\ r\± || r*2j_ (which we call u zz£^) and d±_ _L r*ij_ || r*2j_ 
( u zyy"). We have also measured the orientation-averaged quantity ("ave") denned as 

C a E ve (0;z±; \Ri±\, \R 2 ±\) = J dR 1± J dR 2± C E (9;z ± ;Ri±,R2±), (3-4) 

where J dRi± stands for integration over the orientations of Ri±. The lattice version of this 
equation is easily recovered for even (integer) values of the transverse sizes; in our particular 
case, |rjj_| = 1, we have to use a sort of "smearing procedure, averaging nearby loops as 
depicted in Fig. [2j 

§4. Numerical results and prospects 

In RefsP^ 1 ''^' 1 we have performed a Monte Carlo calculation of the correlation function 
Ql of two Wilson loops for several values of the relative angle, various lengths and different 
configurations in the transverse plane, on a 16 4 hypercubic lattice with periodic boundary 
conditions. The link configurations were generated with the usual Wilson action for SU(3) 
pure-gauge theory, also known in the literature as the quenched approximation of QCD, which 
consists in neglecting dynamical fermion loops by setting the fermion matrix determinant to 
a constant. 

We have measured the correlation functions (Wli>Vl2) and the loop expectation values 
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(Wli) and (W L 2), with W L i =W L {h\\;rir,d) and W L2 = W L {h\\; r 2 r, 0), on 30000 ther- 
malised configurations at = 6/g 2 = 6.0. As it is well known, the lattice spacing a is 
related to the bare coupling constant g (i.e., to (3) through the renormalisation group equa- 
tion. The lattice scale, i.e., the value of a in physical units, is determined from the physical 
value of some relevant (dimensionful) observable like the string tension or the static qq force 
at some fixed distance: in our case one finds that a(/3 = 6.0) ~ 0.1 fm. The choice of /3 = 6.0 
on a 16 4 lattice is made in order to stay within the so-called "scaling window 11 : in this sense 
we are relying in an indirect way on the validity of the relation (I3-3P between Wilson-loop 
correlation functions on the lattice and in the continuum (and therefore we shall use the 
notation Qe/Ce of the continuum in all the figures reporting our lattice data). 

As explained in Section 2, we are interested in the T — > oo limit and so we have to 
somehow perform it on the lattice. In practice, we have to look for a plateau of the correlation 
function plotted against the loop lengths Li and L 2 : in Fig. Owe show the dependence of 
the correlator on the length L\ = L 2 = L of the loops at 9 = 90°. Of course, on a 16 4 lattice 
it is difficult to have a sufficiently long loop while at the same time avoiding finite size effects 
and at best we can push the calculation up to L = 8; nevertheless, a plateau seems to have 
been practically reached at about L = L p \ ~ 6 4- 8. As 9 varies from 90° towards 0° or 180°, 
we expect L p \ to grow. Indeed, L p \ blows up at 0°, 180° due to the relation between the 
correlation function and the static dipole-dipole potential Vdd (see RefsP^l 1 '^ and references 
therein: some preliminary lattice data for Vdd have been obtained in Ref.^): 



0;T;z±;R 1± ,R 2± ) ~ exp 

T— yoo 



-2T Vdd(z±, R\±, R 2 ±) 



1, (4-1) 



from which we expect to diverge at 9 = 0° and, by virtue of the crossing-symmetry 
relations^ also at 9 = 180°. In the following we will consider only 9 ^ 0°, 180°: our 
data show that the correlation function is already quite stable against variations of the loop 
lengths at L\, L 2 ~ 8 (at least for 9 not too close to 0° or 180°) and so we can take the data 
for the largest loops available as a reasonable approximation of Cl, defined as the asymptotic 
value of Ql as L 1; L 2 — > oo. 

We have considered the values d = 0,1,2 for the distance between the centers of the 
loops: as expected, the correlation functions vanish rapidly as d increases, thus making the 
calculation with our simple "brute force" approach very difficult at larger distances. 

From now on we will discuss the issue of the angular dependence of the correlation 
function. As already pointed out in the Introduction, numerical simulations of LGT can 
provide the Euclidean correlation function only for a finite set of ^-values, and so its analytic 
properties cannot be directly attained; nevertheless, they are first-principles calculations that 
give us (inside the errors) the true QCD expectation for this quantity. Approximate analytic 
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Fig. 3. Dependence of Ge on the length L\ = 
L2 = L (in lattice units) of the loops at 
= 90° ford = 0,1, 2. 



Fig. 4. Comparison of the lattice data to the 
SVM prediction (|4-2|) with -K"svm calcu- 
lated according to Ref.® (solid line) and 
to the one-parameter (.Ksvm) best-fit (for 
the u zz£' and a zyy" cases only) with the 
SVM expression fl32J) (dotted line) at d = 
1. 



calculations of this same function then have to be compared with the lattice data, in order to 
test the goodness of the approximations involved. The Euclidean correlation functions we are 
interested in have been evaluated in the Stochastic Vacuum Model (SVM),® in the Instanton 
Liquid Model (ILM) ™™ and using the AdS/CFT correspondence the comparison of our 
data with these analytic calculations is not, generally speaking, fully satisfactory. 
In the SVM® the Wilson-loop correlation function is given by the expression 

C£ YM) {6) = ~ exp (-^K SVM cot 0\ + i exp f~# S VM cot ej - 1, (4-2) 

where -Ksvm is a function of z±, Ri± and R\±_ only, whose precise expression, given in Ref.,® 
we have used to numerically evaluate the correlator (14-21) in the relevant cases. The SVM 
prediction (14-21) agrees with our lattice data in a few cases, at least in the shape and in 
the order of magnitude, but, in general, it is far from being satisfactory: for example, the 
comparison with our data for d = 1 is shown in Fig. HJ More or less the same conclusion is 
reached if one instead performs a one-parameter (ITsvm) best-fit with the given expression: 
the values of the chi-squared per degree of freedom (Xd.o.f.) °f this and the other fits that 
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Table I. Chi-squared per degree of freedom for a best-fit with the indicated function. 



Xd o f 


d = 







d=l 






d = 2 






zzz/ zyy 


ave 


zzz 


zyy 


ave 


zzz 


zyy 


ave 


SVM 


51 




16 


12 




1.5 


2.2 




pert 


53 


34 


16 


13 


13 


1.5 


2.2 


4.5 


ILM 


114 


94 


14 


15 


45 


0.45 


0.35 


1.45 


ILMp 


20 


9.4 


0.54 


0.92 


1.8 


0.13 


0.12 


0.19 


AdS/CFT 


40 




1 


0.63 




0.14 


0.065 





we have performed are listed in Table [H 

We have also tried best-fits with the following simple functional forms: 

c < i eTt \e) = K pcit (cote) 2 , 



C 



(ILM) 



(0) 



sin 9 



(4-3) 
(4-4) 



The first expression ( 14- 3p is exactly what one obtains in leading-order perturbation the- 
ory [ID JED J1D r-Q^ second expression ()4-4p is the one predicted by one-instanton effects in the 
ILM p,[l9} The resu its, shown in Table H are again not satisfactory. In particular, the ILM 
expression seems to be strongly disfavoured at d — 0, while at d — 2 it looks better than the 
SVM and perturbative-like expressions. 

By combining the two previous expressions into the following expression, 



C 



(ILMp) 



(4-5) 



largely improved best-fits have been obtained, as one can see in Table HI The resulting 
best-fit functions in the d = 1 cases are plotted in Fig. [5J 

A well-defined numerical prediction for the prefactor Kjlm of 1/ sin 9 in the ILM has been 
obtained in Ref.P in Table O we compare this prediction with the value obtained with a 
fit to the lattice data with the fitting functions (14-41) and (14- 5p . The ILM prediction turns 
out to be more or less of the correct order of magnitude in the range of distances considered, 
at least around 9 = tt/2, but it does not match properly the lattice data. The agreement 
with the data seems to be quite good at d — 2; however, concerning the dependence on 
the relative distance between the loops, it seems that the ILM overestimates the correlation 
length which sets the scale for the rapid decrease of the correlation function. This is also 
supported by the comparison of the instanton-induced dipole-dipole potential Vm with some 
preliminary numerical results on the latticed 

Finally, we have tried a best-fit with the expression that one obtains using the AdS/CFT 
correspondence, for the M = 4 SYM theory at large iV c , large 't Hooft coupling and large 
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Fig. 5. Comparison of lattice data to best-fits 
with the perturb at ive-like expression (|4-3p 
(solid line), the ILM expression (|4-4p (dot- 
ted line) and the ILMp expression (|4-5p 
(dashed line) at d = 1. 



Fig. 6. Comparison of lattice data to a best- 
fit with the AdS/CFT expression (|4-6p for 
various cases. 



Table II. Value of Kilm X 10 3 for the relevant configurations: ILM prediction (first column), "ILM" 
fit (second column) and "ILMp" fit (third column). 





predicted 


fitted-ILM 


fitted-ILMp 


d 


zzz zyy 


zzz zyy 


zzz zyy 





0.880-1.08 0.880-1.08 


17.1 17.1 


9.86 9.86 


1 


0.827-1.02 0.798-0.984 


7.60 5.79 


4.32 3.39 


2 


0.692-0.853 0.607-0.748 


1.31 1.43 


0.947 1.14 



distances between the loopsP3 



C 



(AdS/CFT) 



(9) = exp 



sin 9 



+ K 2 cot 9 + K 3 cos 9 cot 9 > - 1 



(4-6) 



The results are shown in Table |T] and in Fig. El Taking into account that this is a three- 
parameter best-fit, even this one is not satisfactory: best-fits with QCD-inspired expressions 
with only two parameters, like, e.g., the ILMp expression (14- 5 p [or some appropriate modi- 
fication of the SVM expression (I4-2|) ] give smaller Xd.o.f.- 

As we have said in the Introduction, the main motivation in studying soft high-energy 
scattering is that it can lead to a resolution of the total cross section puzzle, so it is worth 
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discussing what the various models have to say on this point. Using Eqs. (12 -ip . (12-lip and the 
optical theorem, it is easy to see that the SVM and the ILMp expressions give constant cross 
sections at high energy, as in these cases the high-energy limit can be carried over under 
the integral sign, so that the knowledge of the ^-dependence of the correlation function is 
sufficient to completely determine, after the analytic continuation 9 — > —ilog(s/m 2 ), the 
high-energy behaviour of total cross sections. Although the AdS/CFT expression ( 14 -6p is 
not, of course, expected to describe real QCD, it nevertheless shows how a non-trivial high- 
energy behaviour could emerge from a simple analytic dependence on the angle 9. In this 
case, after the analytic continuation into Minkowski space-time, it is not possible to pass 
to the high-energy limit under the integral sign, as the integrand is an oscillating function 
of the energy, and one should carry over the remaining integrals first. The experimentally 
observed universality in the high-energy behaviour suggests that the integration over the 
distance between the loops should be the relevant one: this seems to be the case also in the 
AdS/CFT case, where, combining the knowledge of the various coefficient functions in ( I4-6P 
in the large impact-parameter region^ with the unitarity constraint in the small impact- 
parameter region, a variety of possible high-energy behaviours for the total cross section is 
shown to emerge (including, e.g., a pomeron-\ike behaviour a ~ s 1 / 3 ).^ 
It seems then worth investigating further the dependence of the correlation functions on the 
relative distance between the loops, as well as on the dipole sizes, as they could combine in 
a non-trivial way with the dependence on the relative angle: these and other related issues 
will be addressed in future works. 

As a final and important remark, we note that our data show a clear signal of C- 
odd contributions in dipole-dipole scattering. Because of the crossing-symmetry relations 
(I2-I2p . it is natural to decompose the Euclidean correlation function Ce{9) as a sum of 
a crossing-symmetric function C~^{9) and a crossing-antisymmetric function 0^(9), with 
Cg(9) = | [Ce(9) ±Ce(tt — 9)]^ Upon analytic continuation from the Euclidean to the 
Minkowskian theory and using Eq. (12-121) . one can show that they are related respectively to 
pomeron (i.e., C = +1) and odderon (i.e., C = —1) exchanges in the dipole-dipole scattering 
amplitude. Looking at our lattice data, we notice that there is an asymmetry with respect to 
9 = 7r/2 in the plot of the Euclidean correlation function, for the u zzz" and ll zyy n tranverse 
configurations, against the relative angle. (As regards the orientation-averaged quantity 
C E ve ', defined in Eq. (13-41) . it is trivially crossing-symmetric by virtue of Eqs. (12-121) .) 
In other words, a small but non-zero crossing-antisymmetric component is present in 
our data, thus signalling the presence of odderon contributions to the loop-loop correlation 
functions and in turn to the dipole-dipole scattering amplitudes. Even though these C-odd 
contributions are averaged to zero in meson-meson scattering (at least in our model, as long 
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as the squared meson wave functions satisfy some reasonable symmetry properties in their 
dependence on the dipole orientations and on the longitudinal-momentum fractions), they 
might play a non-trivial role in more general hadron-hadron processes in which baryons and 
antibaryons are also involved. 
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